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Abstract. 

In this talk we present a new method for determining the nonperturbative quark-gluon 
vertex, which constitutes a crucial ingredient for a variety of theoretical and phenomenological 
studies. This new method relies heavily on the exact all-order relation connecting the 
conventional quark-gluon vertex with the corresponding vertex of the background field method, 
which is Abelian-like. The longitudinal part of this latter quantity is fixed using the standard 
gauge technique, whereas the transverse is estimated with the help of the so-called transverse 
Ward identities. This method allows the approximate determination of the nonperturbative 
behavior of all twelve form factors comprising the quark-gluon vertex, for arbitrary values of 
the momenta. Numerical results are presented for the form factors in three special kinematical 
configurations (soft gluon and quark symmetric limit, zero quark momentum), and compared 
with the corresponding lattice data. 


1. Introduction 

One of the major challenges of nonperturbative QCD is to understand the mechanism that 
drives chiral symmetry breaking and the associated generation of constituent quark masses. As is 
well known from a series of previous studies [n El Elia [5], the quantity that is intimately related to 
the underlying dynamics of the chiral symmetry breaking is the quark-gluon vertex. In addition, 
this vertex is of paramount importance in the formalism of the Bethe-Salpeter equations (BSEs), 
which describes the formation of the bound states of the theory 0 El El El EDI El]. 

From the point of view of the perturbative QCD, the quark-gluon vertex has been carefully 
scrutinized at the one-loop level [TD], where results for general kinematic configurations and 
arbitrary gauges. Moreover, at two and three-loop order we have results for some specihc 
kinematics and gauges [laEii. 

The main difficulty in dealing with the quark-gluon vertex lies in the fact that one has to 
determine the behavior of twelve form factors (four “longitudinal” and the eight “transverse”), 
which are functions of three independent kinematic variables. For this reason, the available 
nonperturbative information on this quantity is rather limited. In particular, there are only 
few results obtained from simulations on relatively small lattices da EH EZl Ea EH ED]. In 
the context of the Schwinger-Dyson equations (SDFs) the situation is not that different. The 
behavior of the form factors is governed by a complex system of coupled integral equations, 
which can be solved only after considerable truncations and drastic simplifying assumptions 

[HlEllEIlEaEaEaES]- 


In this talk we will present a novel nonperturbative approach for calculating the form factors 
of the quark-gluon vertex in the Landau gauge m- This task will be accomplished within the 
PT-BFM scheme [281 [29l [30], which is obtained from the combination of the pinch technique 
(PT) [3ll[32l[33l[3lll35l[36| with the background field method (BFM) [37]. 

An intrinsic feature of the PT-BFM formalism is the natural separation of the gluonic 
field into a quantum and a background part, thus increasing the number of possible Green’s 
functions that one may consider. In particular, two types of the quark-gluon vertices make 
their appearance: (i) the conventional quark-gluon vertex (formed by a quantum gluon, quark 
and anti-quark fields), denoted by F; and the background quark-gluon vertex (formed by a 
background gluon, quark and anti-quark fields), denoted by F. A crucial difference between these 
two vertices, lies in the fact that the conventional vertex satisfies the usual STIs, whereas the 
background vertex obeys Abelian-like WIs. The conversion between quantum and background 
vertices is achieved through the so-called background-quantum identities (BQIs) [38ll39], which 
relate F and F through special auxiliary ghost Green’s functions, namely A^^,, and its 
conjugated K^. 

The aim of this work is to express the conventional quark-gluon vertex as a deviation from 
the “Abelian-like” vertex F^. Specifically, our strategy will be the following: first we use the 
“gauge technique” inspired Ansatz gomiiiiiiiis] for the longitudinal part of the Abelian-like 
F^. Then, with the help of the so-called “Transverse Ward Identities” (TWIs) [SI IT6l ITTIBS] . 
we will fix the transverse part of the F^, neglecting the non-local terms present in the TWIs. 
The combination of these two steps generates the so-called minimal Ansatz for F^ [S]. As a 
final step we use the BQIs to convert F^ into F. This conversion requires the computation of 
the special auxiliary functions and the behavior ofA^j^ is well-known from previous 
studies [50], whereas and were first computed at one-loop dressed approximation in 
Ref. [27], using as ingredient the gluon lattice propagator. 

In order to make contact with the results of lattice simulations [151 dsi da dsi dsi [20] , we 
will present the numerical evaluation of the relevant form factors in three special kinematical 
configurations namely (i) soft gluon, (ii) quark symmetric limit and (iii) zero quark momentum. 
As we will see, while a general qualitative agreement with the available lattice data is found, 
quantitative discrepancies still remain. 


2. The two quark-gluon vertices in the PT-BFM formalism 


In the PT-BFM formalism there are two quark-gluon vertices, depending on the nature of the 
incoming gluon. Specifically, the vertex formed by a quantum gluon (Q) entering into a pah 
corresponds to the conventional vertex, F“ (see left vertex of Fig. dj; the corresponding three- 
point function with a background gluon (A) entering represents instead the PT-BFM vertex, 
and will be denoted by F" (see right vertex of Fig. [1]). Ghoosing the flow of the momenta such 
that Pi = q + P2, we then define 


*r“(g,P2, -Pi) = igt‘"T^{q,p2, -Pi)] *r“(g,p2, -Pi) = igt''fii,{q,P2, -pi), (1) 


where the hermitian and traceless generators of the fundamental SU(3) representation are 
given by = A“/2, with A“ the Gell-Mann matrices. 

It is important to stress that F^ and F^ coincide only at tree-level, where one has 
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The essential difference between these two vertices is that F^ obeys the QED-like WI [37| 


q'"'riiiq,P2,-Pi) = S ^{pi)-S \p2), 


( 2 ) 
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Figure 1. The conventional and background quark-gluon vertex with the momenta routing 
used throughout the text. 


instead of the standard STI 

-Pi) = F{q^) [S~^{pi)H{q,p 2 , -pi) - H{-q,pi, -p 2 )S~^{p 2 )\ , (3) 

satisfied by T^. In both expressions, S~^{p) = A{p‘^)'^ — B{p‘^), is the inverse of the full quark 
propagator, defined in terms of the wave function, A{p'^), and the mass function B{p^). In 
addition, F{q^) denotes the ghost dressing function, related to the full ghost propagator by 
D{q^) = F{q^)/q^, whereas the functions H°^ = —gFH and its conjugated = gFH 
correspond to the so-called quark-ghost kernel, and are shown in Fig. [2j 

“ i i'' 
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Figure 2. The ghost kernels H and H appearing in the STI satisfied by the quark vertex T^. 
The composite operators and ijjc^ have the tree-level expressions —gf^ and gF respectively. 


Notice that the quark-ghost kernel admits the general decomposition [T2] 

H{q,P2, -Pi ) = + -^ 1^1 + -^ 2^2 + X 3 a^uPiP 2 ) ( 4 ) 

Fi-q,pi,-p2) = XoI +X2:^i + Xi;^2 + X3afj,i,p^P2, 

where 7i/] and we have introduced the compact notation on the form factors 

Xi = Xi{q^,p^,p\) and Xi = Xi{q^,p\,p^). Notice that at tree-level, one clearly has 
xj^^ = X® = 1, with the remaining form factors vanishing. 






It is important to stress here that a set of identities, called background quantum identities 
(BQIs) [381139| ■ relate the conventional and PT-BFM vertices. The BQI of interest in our case 
reads |30] 


r^(9,P2,-Pi) = 






^v{q,P2, -pi) 


- S ^{pi)K^{q,p2,-pi) - K^,{-q,pi,-p2)S ^(^ 2 ), 


-1/ 


( 5 ) 


where, in the Landau gauge, the special functions and its conjugated are related to the 
quark-ghost kernel H (and its conjugated) in the following way 


H{q,P 2 ,-pi) 

H{-q,Pi,-p2) 


l + q^^K^{q,P2,-pi), 

1 - q^K^{-q,pi,-p2). 




a,fi 



iKl{-q,pi, -P2) 



a, f.1 



( 6 ) 


Figure 3. The auxiliary functions and appearing in the BQI relating the conventional 
quark vertex T^ with the PT-BFM vertex F^. The diagrams on the right represent the one-loop 
dressed approximation of the two functions. 


In what follows we will use the one-loop dressed approximation, in which the propagators are 
fully dressed while vertices are retained at their tree-level (see Fig. [3] again). This approximation 
yields the following expressions 

K^2iq,P2,-pi) = j^S{k+p2)YPi,y{k)A{k^)D{k- q), 

K^{-q,Pu -P 2 ) = ^9^Ca - k)P^,{k)A{k^)D{k - q ), 


( 7 ) 
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Figure 4. The auxiliary function appearing in the BQI relating the conventional quark 
vertex T^ with the PT-BFM vertex T^. The diagrammatic representation of the gluon-ghost 
scattering kernel 



where Ca represents the Casimir eigenvalue of the adjoint representation [Ca = N for 
SU(N)], d = 4 — e is the space-time dimension, and we have introduced the integral measure 
= d'^A;/(27r)‘^, with // the’t Hooft mass. 

The functions G{q‘^) and L{q^) appearing in Eq. ([2]) are the tensorial projections of the special 
two-point function 


= -ig'^CA ^J,{k)D{q - k)Hy^{-q,q - k,k) 

J k 

= 9^iuG{q^) + ( 8 ) 

qz 

Finally, H^y is the ghost-gluon scattering kernel, and A^y{q) is the gluon propagator defined in 
the Landau gauge as 

iAfj,y{q) = —iP^y{q)A{q ), P^y^q) = g^y — q^qy/q ■ (9) 

Interestingly enough, in the Landau gauge the form factors G{q^) and L{q^) are related to 
the ghost dressing function F{q^) by the all-order relation [CTl [HU] 

F-\q^) = l + G{q^)+L{q^). (10) 

The most general tensorial decomposition of the quark-gluon vertex contains four longitudinal 
and eight transverse form factors. Using transverse and longitudinal (T-|-L) basis [121152| . we 
can write 


r''(g,P2,-Pi) = '^^Kq^^PhPi)Li{q,P2,-Pi) + '^^f{q^,php\)Tl'{q,P2,-pi), (11) 

^=1 i=l 

where the longitudinal basis vectors read (remember that t = pi + P 2 ) 




L>i = a>^^ty; 


( 12 ) 



while for the transverse basis vectors we have instead 


Ti" =p^{pi ■ q)-Pi{p2 ■ q)] 

Ty = -^{q-t)L^ - tf^a,yxPiP2'^ 


= Tf/; 

= T^^.xPiPh 

Ti:=r{q-t)-t^-, 

= -yf^a^xPiP^ + P2^i - Pih■ (13) 


In addition to the usual WI ([21) and STI ([3]), specifying the divergence of the quark-gluon 
vertex there exists a set of less familiar identities, called transverse Ward identities (TWIs) 

gaiiaiMiiiziiiEiiiH], which give information on the curl of the vertex, d^Ty — dyV 

In the case of an Abelian gauge theory, one may consider a fermion coupling to a gauge boson 
through a vector vertex and an axial-vector vertex In this case the TWIs read |49] 


qv^u{q,P2, -pi) - qu^^l{q,P2, -pi) = i[S~^{p2)'^^lu - -b 2imT^y{q,p2, -pi) 

+ t^ex^ypTi{q,p2, -pi) + A'^y{q,p2, -pi), 
qu^ti<i^P2, -Pi) - q,y^p{q,P2, -pi) = i[S-^iP2)^py - af^yS-^{pi)] 

+ t^expypTP{q,p2, -pi) + V^y{q,p2, -pi), (14) 

where a^^y = j^a^y, cxpyp is the totally antisymmetric Levi-Civita tensor, while T^y, AJ^y, and 
Vpy represent non-local tensor vertices that appear in this type of identities. 

As was shown in Ref. |49] . it is possible to disentangle the vector and the axial-vector vertices 
appearing in Eq. HUD. Specihcally, for the vector vertex we obtain 

[tf^ei^qp - {q-t)ei]rP{q,p2, -pi) = Pr{i[S~^{P2)'5py - ^pyS~^{pi)] + V^y{q,p2, -pi)}, (15) 


where the tensorial projectors are defined as 


= * = 1 , 2 ; el = t^, = ( 16 ) 

In what follows we will use Eq. (jl5p in conjunction with the WI ([2]) in order to determine 
the complete set of form factors characterizing the vertex E^. After that, we will apply the BQI 
given by Eq. Q to obtain the hnal expression for E^. 


3. Special kinematic configurations 

In this section, we will present the results for three special kinematic configurations: (i) the 
soft gluon limit, (ii) the symmetric limit and (hi) the zero quark momenta configuration. The 
general expressions, which are valid a generic configuration, can be found in Ref. m- 


3.1. Soft-gluon limit 

Let us start with the soft limit, obtained when we take the limit of pi —)• p 2 or similarly 
g —)■ 0. In this limit all the transverse tensor structures (jl3p vanish identically. The vertex is 
therefore purely longitudinal, and after setting pi = P 2 = P, the longitudinal tensorial structures 
reduce {pi = P 2 = p) to 


L^i = 7^; 




= 2p^; 


L>f = 2a>^^py. 


(17) 


The form factors that accompany each one of the above tensorial structures are given by 

Fq^TI = A{1- 2p^K^) - 2BK^, 

Fq^T 2 = 2 A' + 2 A (K 3 + K 4 ) - 2 BK 2 , 

FQ-^rS = -2B' + 2A (K^ + p^K 2 ) - 2 BK 3 , 

r4 = 0, (18) 

where Fg"^ = F“^(0), A = A(p^), B = B{p^), = K-{p‘^), and a prime denotes the derivative 

with respect to p‘^. The K- correspond to the decomposition of the tensorial structure of Eq. © 
in the basis presented in Eq. m, and its detailed derivation is given in Ref. EZj. 

All ingredients that are necessary for computing the K- = K^{p^) are renormalized at 
p = 2.0 GeV. In particular, we use the SU(3) gluon propagator, A((jr), obtained by the lattice 
simulation of Ref. m, the solution of the SDE for ghost dressing function, F(q), and the 
auxiliary functions 1 + G{q) and L{q). All these quantities where computed using a{p) = 0.45. 
In addition, the behavior of the functions A{p) and B{p) were obtained by solving the quark 
gap equation for a current mass niQ = 115 MeV. 

In Fig. [5] we plot the functions Ki in the soft gluon limit. With the Ki at hand, the next step 
is to determine the vertex form factors of Eq. (USD. 



Figure 5. (color online). The auxiliary functions Ki evaluated in the soft gluon limit. 
Specifically, in Figs. Eland [7] we plot the form factors 

Ai(p) = rf(pE); A2 (p) = ^E^Pe); A3(p) = -^TKpe), (19) 

and compare them with the lattice data of m, obtaining rather reasonable agreement. 


3.2. Symmetric limit 

The symmetric limit is dehned taking pi ——P 2 - In this limit, only one longitudinal basis 
tensor m and two transverse tensors m survive. Specifically we have 


'^3 = 4 - p '"#); 


= r-, 


= -25^V- 


( 20 ) 














Figure 6. (color online). The soft gluon form factors Ai (left) and pX^ (right). Lattice data in 
this and all the following plots are taken from m- 



Figure 7. (color online). The form factor A 2 and the corresponding lattice data. 


However, the lattice simulations can not determine separately the form factors that 
accompanies the tensorial structures above described; only combinations of the type 

g2p{n +P^^l) = + ^(1 + ‘2p^Kl) - 2B{Ki+p^Kl), 

g2pTl = 2B' + 2A (iLf + _ 2BK^. (21) 

can be extracted. In the above equation the compact notation A = A{p'^), B = B{p^) and 
= 1 + G{Ap^) has been introduced. 

Using the same ingredients described in the soft gluon case, we have computed the 
corresponding Ki for the symmetric configuration, which are presented in Fig. [51 In this limit, 
we clearly see that we have a divergent K 2 and a finite iL 4 . 

The next step is to compare our numerical results for the Euclidean version of the form factors 
combinations of Eq. (1211) with the lattice data of Ref. m- This comparison is shown in Eig. [9l 
where we have defined 

Kip) = ^i'^ipK -pl'^riPB); nip) = ^rp(pB), 


(22) 





















































Figure 8. (color online). The auxiliary functions Ki evaluated in the symmetric gluon limit. 



Figure 9. (color online). The symmetric limit form factors and T 5 compared with the 
corresponding lattice data. The grey curves are obtained through simple rescaling of the blue 
ones. 


Clearly, we see that the overall shape of the both form factors are correctly reproduced; 
however, the overlap with the lattice data is only obtained if we rescale our results by (different) 
multiplicative factors, giving rise to the gray curves. 

In Figs. [To] and m we plot, for completeness the three non-zero form factors separately. 
Notice that both T]" and Tg are finite, whereas Tg is divergent. 

3.3. Zero quark momentum 

The next limit that we will present here is the so-called zero quark configuration, where we 
set to zero the quark momentum p 2 , which leads to q = pi = p. A crucial difference of this case 
compared with the previous one is that the Kf’^{p^,t),p^) and K^' (p^,p^, 0 ) do not coincide 
anymore, and need to be evaluated separately. More specifically, the non-zero form factors are 
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Figure 10. (color online). The symmetric form factors T]" (left) and 
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Figure 11. (color online). The divergent form factor Tg in the symmetric limit. 


expressed as 

F-^ri = A{l+p^K^) - BK^ - BoKi, 

= - ^^3 - BqKI 

A{l+p^Ki)-BK^ , 

gVl = -^{B - Bo) - A{K{ +p^Kl) - BKl - BoKI, (23) 

with the usual definitions ^4 = ^(p^), B = B{p‘^), as well as Bq = B{0). 

In Fig. [12] we show the auxiliary functions K- (left) and (right) evaluated in the zero 
quark momentum configuration. Notice that all K- and K^ are IR finite except for K 2 and K 2 
which are divergent. 

In Fig. [13] we show the Euclidean version of the numerical results for the form factors 
appearing in Eq. (I23|) . In particular, we see the appearance of a negative divergence in the 
form factor F™. 


grl = -A{Ki + Kl) - BKl - BoKo + 












Figure 12. (color online). The auxiliary functions (left) and (right) evaluated in the 
zero quark momentum configuration. 
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Figure 13. (color online). The form factors Tg® and T™ evaluated in the zero quark 
momentum configuration. 
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Figure 14. (color online). The form factors (left) and rs (right) in the quark zero momentum 
configuration. 




Although there is no lattice results for this kinematic limit, the combination that could be 
measured on the lattice would be of the type 

g(r[ +/r^) = A{l-P^K^) - B{K^+p^Kl) - b^{k\+p^kI), 

gvl = Bo) - A{K[ + p^Kl) - BKl - Bol^,. (24) 

In this case, we present our “prediction” for these combinations in Fig. [TTl 
4. Conclusions 

We have presented a new methodology for determining the longitudinal and transverse form 
factors of the nonperturbative quark-gluon vertex within the PT-BFM scheme. This scheme 
allows us to take the full advantage of the rich amount of information originating from the 
fundamental underlying symmetries, which are encoded in a set of crucial identities such as 
WIs, STIs and BQIs. 

The key observation in this analysis is the connection between the two distinct quark-gluon 
vertices, T^ and T^, appearing in the PT-BFM scheme. Using the WIs and the TWIs satisfied 
by P^, we first determine the form factors that describes the behavior of this Abelian-like type 
of vertex. Then, with the help of the BQI that connects both vertices, we obtain the final 
expression for the conventional P^. 

We have shown that the BQI is expressed in terms of the auxiliary three-point functions 
which were calculated in the one-loop dressed approximation. Already, at this level of 
approximation, we have obtained nontrivial information for all form-factors. In addition, we 
have noticed that the contributions originating from the K- and are in general sizeable, and 
therefore the and contribute significantly in obtaining results similar to those found in 
lattice simulations. 

For the determination of the K- and we have used as external ingredient the full gluon 
propagator, A(g^), obtained in lattice simulations. The remaining necessary ingredients, namely 
the ghost dressing function F{q^) and the quark functions A{p^) and B{p‘^), were obtained 
solving numerically their corresponding SDEs. 

For the purpose of this talk we have applied our formalism to three particular kinematic 
limits known as (i) “soft gluon”, (ii) “quark symmetric” and (ii) “zero quark momentum” 








configurations, which give rise to considerable technical simplifications, especially in the 
calculation of the and K^. Evidently, the numerical analysis presented here may be extended 
to arbitrary kinematic configurations, furnishing valuable information on such a fundamental 
quantity as the quark-gluon vertex, which constitutes a crucial ingredient for a variety of 
theoretical and phenomenological studies. 
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